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ABSTRACT 

Tests of big bang nucleosynthesis and early universe cosmology require preci- 
sion measurements for helium abundance determinations. However, efforts to de- 
termine the primordial helium abundance via observations of metal poor H II re- 
gions have been limited by significant uncertainties (compared with the value in- 
ferred from BBN theory using the CMB determined value of the baryon density). 
This work builds upon previous work by providing an updated and extended 
program in evaluating these uncertainties. Procedural consistency is achieved 
by integrating the hydrogen based reddening correction with the helium based 
abundance calculation, i.e., all physical parameters are solved for simultaneously. 
We include new atomic data for helium recombination and collisional emission 
based upon recent work by Porter et al. and wavelength dependent corrections to 
underlying absorption are investigated. The set of physical parameters has been 
expanded here to include the effects of neutral hydrogen collisional emission. It 
is noted that H7 and H<5 allow better isolation of the collisional effects from the 
reddening. Because of a degeneracy between the solutions for density and tem- 
perature, the precision of the helium abundance determinations is limited. Also, 
at lower temperatures (T < 13,000 K) the neutral hydrogen fraction is poorly 
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constrained resulting in a larger uncertainty in the helium abundances. Thus the 
derived errors on the helium abundances for individual objects are larger than 
those typical of previous studies. Seven previously analyzed, "high quality" H II 
regions are utilized for the primordial helium abundance determination. The 
updated emissivities and neutral hydrogen correction generally raise the abun- 
dance. From a regression to zero metallicity, we find Y p as 0.2561 ± 0.0108, 
in broad agreement with the WMAP result. Tests with synthetic data show a 
potential for distinct improvement, via removal of underlying absorption, using 
higher resolution spectra. A small bias in the abundance determination can be 
reduced significantly and the calculated helium abundance error can be reduced 
by ~ 25%. 

Subject headings: HII Regions: abundances — galaxies: abundances — cosmol- 
ogy: early universe 



1. Introduction 



Big bang nucleosy nthesis (BBN) provides us with a unique window to the Uni verse at 



redshifts of order 10 10 (iWalker et al.l 



1991 



Olive et al. 



2000 



Fields fc Sarkar 



2008). With 



the WMAP determination of the baryon density of the universe feomatsu et al.l 120091 ) 



BBN has become a zero-parameter theory ( jCyburt et al. 



2002 



20031 ). leadin g to relatively 



preci s e predictions of t he lig h t element abundances of D, 3 He, 4 He, and 7 Li (ICvburt et al. 



2001 



Coc et al. 



Serpico et al 



2002 



2004 



2004: 



Durt 



Cyburt et al 



2004 



Descouvemont et al. 



12004 : 



ICuoco et al. 



2004 



20081 ). In addition to lending itself to scrutiny by com- 



paring the predictions of BBN with observations, the resulting abundances are of significant 
scientific value in constraining the cont e nt and interactions of the universe during BBN 



( IMalaney &: Mathews 



1993 



Sarkar 



1996 



Cyburt et al.ll2005l ). These constraints, however, 



generally require accurate determinations of the light element abundances, particularly in the 
case of 4 He where better than 1% determinations are necessary. The difficulties in obtaining 
this precision for 4 He is the focus of this paper. 

The concordance between BBN and observations of the CMB temperature anisotropies 
is predominantly due to the agreement between the BBN predicted primordial value of D /H 
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compared with determinations made from observations of high redshift quasar absorption 
systems. The cosmic microwave background (CMB) power spectrum provides a high pre- 
cision measurement of the baryon density or equivalently the baryon-to-photon ratio, rj, 
(primarily through the first and second acoustic peak heights). Using the 5- year W MAP 
value of tt B h 2 = 0.02273 ± 0.00062 or n = (6.23 ±017) x 10" 10 komatsu et alJbopgh . one 
expects D/H = (2.49 ± 0.17) x 10~ 5 ' (ICvburt et al.ll2008) wh ich should be compared with 
(2.82 ±0.21) x 10~ 5 f Kirkman et all 1200.4 I pettini et al.lboOS, and reference s therein). On 



2008 



Hosford et al. 



the o ther hand, 7 Li is predicted to be substantially higher ( ICvburt et al. 
20091 ) than determinations of 7 Li abundances in the outer layers of halo dwarf stars. This 



places extra emphasis on the determination of the primordial helium abundance, Y p . As a 
result of the very weak logarithmic relationship between r\ and Y p , the measured error on 
Y p must be very small (< 1%) to achieve a reasonable constraint on the determination of 
r\. For comparison, the WMAP value of r\ enables one to predict the value of Y p to within 
0.08%. 

Being relatively chemically unevolved, low metallicity H II regions in dwarf galaxies can 
be used to provide a measure of Y p . Since metallicity only increases with age, these galaxies 
can be used to fit the helium abundance versus metallicity and therefore allows one to extrap- 



olate back to very low values of O/t 



, used here as a surrogate for meta 



licity, corresponding to 

the primordial helium abundance (IPeimbert fc Torres-Peimbertlll974l ). However, after more 



than three decades of work in determining the abundance using this method, the measure- 
ments frequently have not agreed with each other and until recently, post- WMAP, have not 
been in accordance with the WMAP value (see Figure [J below) , Y p = 0.2486 ±0.0002, which 
is based on the BBN calculation of ICvburt et al.l (120081 ) assuming the WMAP determined 
value of rj. 

The quoted errors in the extragalactic determination of Y p have been historically very 
small in comparison to the range in the results (again see Figured]). The early determi- 
nations of Y p focused on statistical errors in the measured helium line fluxes and, as a 
result, underestimated the error. Ideally, the target quantity, the helium to hydrogen ratio 
(y + = „(ff E //? ) in extragalactic H II regions, would directly follow from an observed helium 
to hydrogen emission line ratio. However, as a result of interstellar reddening, underlying 
stellar absorption, optical depth, and collisional corrections, the actual determination of y + 
is complicated, and the true uncertainty can be large. Interstellar reddening occurs as the 
photons are scattered by dust on their journey. Being a wavelength dependent process, each 
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line is affected differently. The stellar continuum juxtaposes absorption features under neb- 
ular emission lines. The H II region itself can absorb some of the emitted photons as well. 
Finally, the total emission fluxes are due to both recombination and collisional excitation. 
In the case of the hydrogen emission, only the small fraction of recently recombined neutral 
atoms can contribute to collisional emission. Unfortunately, these confounding processes can 
not be directly measured. In addition, the temperature, which determines the physical rates, 
cannot be determined directly and independently unless it is assumed to be perfectly uniform 
throughout the H II region. Our ignorance of these necessary physical model parameters will 
necessarily inhibit the determination of y + and therefore correspond to potentially large in- 
creases in the error on that determination. The error on Y p also varies widely due to sample 
sizes. The error on the intercept, Y p , will shrink with increasing number of regression points; 
however, poor quality spectra do not allow for a robust estimation of the systematic effects 
listed above. Therefore, this work focuses on a small number of "high quality" extragalactic 
objects. 



Previous work by 



Olive fc Skillman in 



"self-consistent" methods of 



Paper I (120011) andlPaoer III ( 120041 ) . following the 



Izotov. Thuan. fc Lipovetskv (1994) in determin i ng th e electron 



density (used for the collisional correction) and iPeimbert. Peimbert. fc Ruizl (120001 ) in deter- 
mining the electron temperature (used for the helium emissivities) , addressed some of the 
systematic effects by solving for physical conditions via a minimization of the differences 
between the values of y + calculated from each of the helium lines. The y 2 minimization 
is accomplished by a down-hill simplex algorithm. Monte Carlo simulations over the input 
fluxes were used to measure the statistical variance with the ultimate effect of further in- 
creasing the total uncertainty. Details of the determination of the model parameters and the 
impact of the Monte Carlo process on these parameters will be discussed further in section 
|2] below. 

The purpose of the current work is to explore further refineme nts to the approach of Pa- 
pers I & II. First, there now exist improve ments to the atomic data (IPorter. Ferland. fc MacAdam 



2007, PFM) that can be incorporated (e.g.. 



Izotov et al 



2007 



Peimbert. Luridiana. fc Peimbert 



20071 . hereafter ITS07 & PLP07 respectively). Second, previously, the flux ratios of the hy- 



drogen lines were used to solve for reddening and underlying hydrogen absorption, while the 
flux ratios of the helium lines were used to solve for density, temperature, underlying helium 
absorption, helium self-absorption, and the helium abundance. There is no need to solve for 
these two sets of parameters separately, and doing so artificially suppresses any potential 
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Fig. 1. — Recent historical progression of the primordial helium mass fraction over the pre- 



vious d ecade and compar i son to the WMAP result (wit 



sion b y IPaeel et al. 
1997 



1992 



Izotov. Thuan. fc 
Izotov. Thuan. fc Lipovetskv Il997l~ 



z 



iipovetskv 



i the values given in date progres- 



Izotov fc Thuan 



19941 lOlive. Skillman. fc Steigman 



1998b 



Izotov et al 



1999: 



Peimbert. Peimbert. fc 




Luridiana et al 
I; ISpergel et al 



Komatsu et al 



2009 



Rui: 



2003 



2000; 



2006; 



Thuan fc Izoto 



Spergel et a 



Izotov et al. 



2003; 



Peimbert 



Izotov fc Thuan 



Peimbert. fc Luridiana 



2004; 



Olive fe S 



ki 



2007; iPeimbert. Luridiana. fc Peimbert 



lman 



2007 



labeled as PSTE92, ITL94, OSS97, ITL97, IT98, ICFGGT99, PPR00, 
TI02, PPL02, LPPC03, WMAP03, IT04, OS04, WMAP06, ITS07, PLP07, WMAP08, and 
AOS09, for this work, respectively) . A general increase in the primordial abundance is 
apparent. Note that, historically, the error bars are typically small relative to the differences 
between studies and with WMAP in particular. 
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degeneracies between them. Thus, here we explore the effects of combining two separate 
"self-consistent" approaches into one. Third, the previous work assumed identical values of 
equivalent width of underlying absorption for all H I or He I emission lines, regardless of 
wavelength, and here we investigate the effects of a wavelength dependence in the underlying 

correcting for the collisional 



stellar absorption (e.g., ITS07). Fou rth, we study the effects o 
emission from neutral hyd : 



Stasihska fc Izotov 



2001 



rogen (e.g., 



Luridiana et al.l 



Davidson fc Kinmanll985 



Skillman fc Kennicutt 



1993|; 



20031 ) . The second change, the integration of the 



hydrogen and helium, is a calculational procedure modification, while the other three are 
improvements in the modeling of the physical processes that impact the observed hydrogen 
and helium emission lines. Detailing the incorporation of these four effects and investigation 
of their impact is the primary aim of this paper and is discussed in sections [2]l6j In sec- 
tion [7J we discuss the importance of ionization correction factors. Section [8] presents a new 
determination of the p rimordia 



observations studied in 



Paper II 



helium abundance resulting from the reprocessing of the 
with comparison to those results, and the results reported 
in ITS07 and PLP07. Section [9] details the effect of the enhancements, remaining errors, and 
possibilities for future improvements. 



2. Determining y H 



As described in 



Paper J . six measured helium emission line fluxes (A3889, 4026, 4471, 



5876, 6678, and 7065) relative to H/3 ( p/hp) ) anc ^ their equivalent widths (W(X)) are used to 
determine the helium abundance. Corrections are made for reddening, underlying absorp- 
tion, collisional enhancement, and radiative transfer: 



F(X) E{HP) 



F(H/3) E(X) w ffgy* /r(A)l + g(A) 



10 /(A)C(ff/3) 



where gfffi is the emissivity ratio of H/3 to the helium line. In the approach of papers 
I and II, the reddening relative to H/3 (C(H/3)) and underlying hydrogen absorption (a#) 
are determined separately by Monte Carlo over the three hydrogen line ratios j™, and 
|f|)ll] The optical depth function f T , and collisional to recombination emission ratio g are 



1 T he reddening coefficients, f(A), are caicuiated from the extinction fits of ICardelli. Clayton, fc Mathis 
(|l989h . 
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both temperature and density dependent (the emissivities are also temperature dependent). 
Therefore the physical parameters needed for the corrections are electron density (n e ), under- 
lying helium absorption (a# e , corrected for in terms of equivalent width), optical depth (r), 
and temperature (T). The parameters are determined simultaneously with the abundance 
by minimizing % 2 , defined as the weighted difference between each helium line's abundance 
and the average. This minimization is performed for each realization in a Monte Carlo over 
the input helium line fluxes and equivalent widths. 

An important result of the Monte Carlo analysis was the discovery of degeneracies in th e 



sense that the fit parameters were not always independently constrained (IPaper IhlPaper It ). 
In particular, the temperature and density exhibited a trade-off whereby increasing the 
temperature while decreasing the density would leave the \ 2 relatively unchanged. The 
reason for this can be seen in the form of the emissivities (see $3]). Four of the lines show 
similar exponential dependencies on temperature and density. As a result, the abundance 
from these lines shifts correspondingly for changes in temperature and density. 

To demonstrate the importance of the degeneracy between temperature and density, 
characteristic model parameters and a helium abundance were chosen and used to generate 
synthetic flux data, \ 2 was then calculated keeping the abundance and all of the traditional 
fit parameters fixed except for the density and temperature. The impact of trading the 
density against the temperature can be clearly seen in Figure [2] which illustrates a well 
constrained quotient of density and temperature but a A% 2 = 2.3 boundary spanning a 
temperature range of 5500 K, a density range of 275 cm -3 , and an abundance variation of 
10%. The shallow, extended minimum results in a large range of best fit parameters for 
the density and temperature upon the Monte Carlo perturbation, and the abundance gains 
a much larger uncertainty (see Figure [3] below). Only by solving for these two parameters 
simultaneously can this degeneracy be discovered. 

The inclusion of wavelength dependent absorption, now an(Hf3) and a# e (A), and neutral 
hydrogen collisional emission, Q(Hf3), will modify the helium to hydrogen ratio as, 

F(X) Em gjjggg 1 1 + 0) 
V F(Hf$) E{\) w V$g*M f T (X) 1 + §(A) • U 

The wavelength dependent absorption is included in order to model variations in the amount 
of stellar absorption and is represented by a relative coefficient for each line. The neutral 
hydrogen collisional emission is included to account for the small population of non-ionized 
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Chi-Squared with y + 




Electron Density [cm -3 ] 

Fig. 2. — A plot of the derived helium abundance as a function of the temperature and 
density for a synthetically generated spectrum. Impressed on this diagram is a contour plot 
of x 2 versus density and temperature for the same spectrum. The synthetic model uses 
n e = 100 cm- 3 , a He = 1.0 A, r = 1.0, T = 18,000 K, and y + = 0.08. The extension of 
the x 2 — 2-3 contour with a strong negative correlation highlights the degeneracy between 
density and temperature (note that for synthetic data Xmin = 0-0)- That the x 2 and abun- 
dance contours are nearly perpendicular demonstrates the impact of the degeneracy on the 
abundance determination (±5%). 
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Fig. 3. — Monte Carlo plot of 1000 solutions based upon synthetic data taken from lPaper 1 



The synthetic model is for n e = 100 cm 3 , a# e = 0.1 A, t = 0.1, and T = 18,000 K. The 
solid circle with error bars marks the original (direct) solution; while the solid square with 
error bars marks the average of the Monte Carlo solutions. Upon performing the Monte 
Carlo, the large range in density solutions (with a corresponding large range in temperature 
solutions), gives a much larger density uncertainty, resulting in a marked increase in the 
abundance uncertainty (1% to 3%). 
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hydrogen which can therefore undergo collisional emission. This process is composed of the 
temperature dependent collisional excitation and the cascade of downward transitions and 
requires the introduc tion of a new parameter, the neutral to ionized hydrogen density frac- 



tion, £. In iPaper III , this was assumed to be negligible. Primarily, the neutral hydrogen 



collisional emission corrects the hydrogen lines but through H/3 shows up in the abundance. 
The integration of the hydrogen reddening and helium abundance calculation is in princi- 
ple accomplished by simply summing their \ 2 contributions and minimizing over the now 
expanded set of physical parameters (n e , a He , r, T, C(H/3), a#, simultaneously. The 
combined \ 2 is given by, 

(X fi (A)-X T (A)) 2 , _(^(A)-y) 2 



X 



E 



where, 



X R (X) 



F(X) 



W(X)+a H W 
W(X) 



y 



W(H/3) 

v y + (A) 



(3) 



(4) 



(5) 



Xfl represents the corrected hydrogen flux, the theoretical hydrogen emission, y the aver- 
age helium abundance, and a is calculated from the measured statistical error and equivalent 
width error. 



In our analysis below, we use the same dataset of H II regions as in IPaper III for calcu- 
lating y + and extracting Y p . Those galaxies were chosen from the llzotov &: Thuanl (ll998bl . 
IT98) sample with screening to minimize uncertainties in underlying absorption, to select 
only low O /H systems, and to avoi d contamination by Galactic Na I a bsorption. In addition 



to thos e seven ta rgets, NGC 346A flPeimbert. Peimbert. fc Ruizll2000l ). which was also exam- 
ined in lPaper III , and I Zw 18 (ITS07) are analyzed. This allows for a comparison with each of 



the five targets used in PLP07. All nine except NGC 346 are in cluded in the recent large sam- 
ple of ITS07. The spectrum for Mrk 1 93 ori ginally reported in llzotov. Thuan. fc Lipovetsky 



(119941 ) and used in 
analyses are reported in 
Mrk 193. 



Olive Sz Skillman 



Izotov et al. 



(2004T) has been reanalyzed and the results of the new 
( 20071 ). Here we use the newly analyzed spectrum for 



Ea ch object is analyzed with the new corrections applied successively to the procedure 



used in 



Paper 11 The analysis is conducted with the parameter space restricted to physically 
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meaningful values, and Monte Carlo average values and standard deviations are returned for 
the model parameters and abundance. Reflecting a straightforward improvement in the 
atomic data and modeling, first the updated emissivities are used in conjunction with sepa- 
rate hydrogen reddening and helium abundance calculations. Next, the hydrogen and helium 
calculations are combined into one integrated calculation. The third change is to introduce 
the slight wavelength dependence in the underlying absorption. Finally, the correction for 
neutral hydrogen collisional emission is incorporated. In the following sections, each of these 
enhancements is first discussed theoretically and then applied to the dataset. 



3. Updated Emissivities 

PFM have recently provided improved helium emissivity fits, as a function of temper- 
ature, and collisional contributio ns, as a function of temperature and density, based upon 
the most recent atomic data (see 



Porter et al. 



20051 ). They rep ort broad agreemen t between 



these emissivities and their plasma simulation code CLOUDY (IFerland et al.lll998l ) of better 
than 0.03% within the temperature range 5000 < T < 25, 000K. However, the abundance 
determination also requires the H/3 emissivity. Though not in the scope of the published 
paper, R. L. Porter (private communication), provided us with the H/3 emissivity in terms 
of the published paper's parameterization, 

fi 5fifiQ x 10 5 

E(H(5) = [-2.6584 x 10 5 - 1420.9 (InT) 2 + 35, 546 InT + — — ] T -1 . (6) 

InT 

In implementing these new emissivities the following difficulty was encountered during 
X 2 minimization. As parameterized by PFM, the density dependence of the emissivities 
weakens with increasing density, ultimately becoming negligible. While this occurs at den- 
sities which are unreasonably high for our objects, the addition of the temperature-density 
degeneracy can and does lead to densities that increase nearly without bound. In prac- 
tice, the Monte Carlo realizations were often found to tend to higher densities that were 
unconstrained by the form of the emissivities (in a region beyond their applicability). The 
resulting \ 2 for these solutions decreased only marginally with the large increase in density; 
so the result is effectively an equivalent solution but which is physically meaningless. Upon 
calculating the Monte Carlo average, even a small number (< 1%) bias the returned values, 
especially the density, to an unacceptable degree. The density dependence in the optical 
depth function disfavors ever increasing density but only if the optical depth itself does 
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not approach zero. Refitting t he em issivities, with the collisional correction, to the form of 
Benjamin. Skillman. fc Smita ( 119991 . BSS) avoids the pitfall while sacrificing less than 1% 



accuracy in the region of interest, 12, 000 < T < 20, 000K and 1 < n e < 300 cm -3 , and 
compares favorably with the BSS fit accuracy] 2 ) Therefore, within the physically applica- 
ble region for our objects, the emissivity values themselves are nearly identical to those of 
PFM. At large densities however, the exponential dependence on the density relegates ever 
increasing densities to ever increasing x 2 ■ 111 effect, the re-parameterized fit protects the x 2 
minimization from unphysical solutions while retaining the num erical acc uracy in the phys- 



ical region. Figures H] and [5] below compare the BSS fits used in lPaper III with the new fits. 
The re-parameterized helium emissivities, including the collisional correction and scaled to 
H/3, are, 



A 88q = o.8779T- a12 ^ 00041 ^ 



E(X) l + g(A) 

'3889 

p m2& — 4 2337 n0 - 085 ~°- 00012ra ' 
_p 4471 = 2 (^lT ' 121-0 ' 00020 "' 
F 5876 = o.754T a212 - - 00051 - 
F 6678 = 2 .639T a244 - - 00054 "< 



F 7065 = 5.903T-°- 519 /(1.462 - (0.127 - .00076n e + 0.000000255n 2 )T) 



3.1. Hydrogen Emission 



In addition t o the updated helium emis sivities, the theoretical hydrogen emission line 



ratios taken from 



Hummer fc Storey! ( 119871 ). were refit to incorporate the reported (very 



weak) density dependence. In principle, this is included for consistency, but became of 
particular interest in physically describing the impact of higher density so as to disfavor 
non-physical solutions. The previous logarithmic form in temperature was retained and 



2 This is different from the approach taken in ITS07 where the BSS emissivities are divided by a linear 
temperature dependent fit to the ratio of the BSS and PFM emissivities. 
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Fig. 4. — Comparison of the PFM emissivities, e(\) > ^° those °f BSS (at n e = 100 cm 3 ). 
The BSS fits are the dashed lines while the PFM fits are solid. The progression is, left 
to right, top to bottom, by wavelength: A3889, 4026, 4471, 5876, 6678, 7065. The PFM 
emissivities plotted here are the refit equations of equation EJ on this plot, the equations 
reported in PFM are within the line thickness of the refit equations. 
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Fig. 5. — The PFM emissivities, , plotted relative to those of BSS (at n e = 100 cm -3 ). 
The newer fits deviate from the old by only several percent, but the relative shifts are clearly 
not the same for all six lines. The three strongest lines, A4471, 5876, and 6678, show similar 
behavior, but 4026 is opposite, and 3889 and 7065 cross the zero deviation line. 
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employed for the density fit. The expanded equations are as follows (with Table [1] listing the 
coefficients), 

X T = J2c tJ (\og(T A )y(\og(D)y (9) 
Table 1: Coefficients for the hydrogen emissivities, cy 
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m 
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2.09E-9 


2.56E-6 
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0.0207 


-1.94E-6 


-6.51E-6 




2 


-0.00770 


4.58E-5 


3.48E-6 



3.2. Tracking the Effects 

A comparison of the resulting He abundances using BSS and PFM is shown in Figure [6j 
The updated emissivities, including the collisional correction, in general raise the calculated 
abundance of helium. However, because the changes in different lines are qualitatively dif- 
ferent (again see Figures @] & [5]) , the update does not categorically raise the abundance. For 
A4471, 5876, and 6678, the H/3 to He emissivity is increased rather uniformly (2-3%) across 
the relevant temperature range in comparison with those of BSS. All three lines exhibit a 
gradual increase with temperature. The behavior of A4026 is similar except that the new 
emissivity is shifted to lower values. A3889 is more complicated in that the emissivity is 
reduced at low temperatures but is enhanced at T > 16,000 K. A7065 is similar, though the 
enhancement begins at 13,000 K. Broadly, A3889, A7065, and A4471 track the BSS emis- 
sivities more closely (within ~ 2%) than the other three lines. Therefore, upon repeating 
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the 



Paper III analysis with the new emissivities, these lines will favor returning a similar 



abundance. The three strongest lines with the smallest relative errors will produce a larger 
abundance, thus raising the average abundance. To decrease the discrepancy between the 
A3889 and 7065 and A4471, 5876, and 6678 abundances, the temperature will decrease, thus 
lowering the average abundance slightly. The weak A4026 is shifted opposite to the three 
strongest lines, but still can have a significant impact as described below. 

As anticipated, in most cases, the updated emissivities of PFM raise the abundance (see 
Figure [6]). This was observed by PLP07 and ITS07. Though, on average, the abundance in- 
creases by ~1%, three of the objects actually exhibit slightly decreased helium abundances. 
The decreased abundances of SBS 940+544N and SBS 1159+545 can be explained by dis- 
crepant high A40 26 fluxes . This has been demonstrated using synthetic data, and using the 
parameters from lPaper III for these two objects, we estimate that A4026 fluxes may be high 
by 15% and 10% respectively. Note that the typical uncertainties in the equivalent widths 
of the A4026 lines and the corrections for underlying absorption are both comparable in 
magnitude to this suspected discrepancy. 

Overall, in comparison to BSS, the PFM emissivities favor a compromise whereby the 
temperature decreases and the average abundance increases. However, as we have seen, this 
effect is not uniform and an overly strong A4026 ca n cause th e abundance to decrease. Thus, 



the regression performed on the seven objectsj of iPaper IJ yields a slightly lower intercept 
as indicated in Figure [6j It is also worth noting that for all nine objects the temperature 
decreases as expected. Since this temperature decrease blunts the abundance increase, cal- 
culations based on temperatures not determined simultaneously with the abundance would 
exhibit a more pronounced abundance increase upon replacing the BSS emissivities with 
those of PFM. Qualitatively very similar to BSS, the new emissivities do not dramatically 
shift the abundances or their pattern. The more accurate atomic data is a distinct improve- 
ment but is not the primary limitation in determining the abundance. 



3 I Zw 18 and NGC 346 d o not pas s the three cuts based on equivalent width, oxygen abundance, and 
radial velocity as discussed in IPaper I] . 
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Fig. 6. — Abundance comparison for the target objects analyzed as in iPaper III with the 
BSS and PFM emissivities. The red thinner lines are for BSS with PFM given in thicker 
black bars. Th e quoted p rimordial helium abundance, Y p , is based on a regression of the 
seven objects of Paper ij (i.e., I Zw 18 and NGC 346 have not been used in the regression). 
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4. Combined H and He 

In Papers I & II, the spectra are analyzed in a two step approach. First the hydrogen 
lines are used in conjunction with the temperature as determined by [O III] emission to 
determine the reddening parameter and underlying hydrogen absorption. These two correc- 
tions are then applied to the helium to H/3 input flux ratios. Subsequently, these corrected 
lines are used to determine the density, helium absorption, optical depth, and temperature 
simultaneously with the abundance. This separation is natural as the dominant hydrogen 
line correction is for reddening. This allows the helium spectra to be corrected for interstellar 
effects before resolving the competing nebular emission and absorption effects on the helium 
lines. However, the independence of these calculations is artificially imposed and could im- 
pair the abundance calculation. Since reddening impacts both the hydrogen and helium 
fluxes (and in the same way), it should be a minimization parameter for both. Furthermore, 
even though the hydrogen emissivity ratios are only weakly temperature dependent, the es- 
timated electron temperature is as physically meaningful for the hydrogen emission as for 
helium. 

In order that the calculations be truly self-consistent, the hydrogen and helium based 
calculations are combined into a single minimization. This has the advantage of being sensi- 
tive to all degeneracies. Conceptually the implementation of combining the two calculations 
is straightforward. Nine line ratios are now input, and all seven physical parameters (n e , 
&He, t, T, C(H/3), &H, £) are determined simultaneously with the abundance via a single 
minimization. Equation |3] explicitly gives this combined \ 2 'i though, for this section the neu- 
tral hydrogen collisional correction is not yet included. The previously observed degeneracies 
provide a further impetus for integrating the hydrogen line based reddening calculation with 
the actual abundance calculation. The larger parameter space allows for new degeneracies 
but the additional calculational effects could break them. 

An important detail of this process is that the deblending of He(A3889) and H8 can now 
be done consistently based upon the solution for reddening and temperature. Previously, in 
the sequential analysis, He(A3889) is deblended using a the hydrogen reddening parameters 
and a fixed temperature. As a result, even though the helium lines are used to determine 
the temperature, the He(A3889) flux could not be adjusted to account for this temperature 
change. This effect is very small since the H8/H/3 emissi yity is only very weakly t emperature 



dependent ( = 0.104 T^- 046 , from a fit to the data of lrlummer &: Storeylll987 ). However, 
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integration allows the blended input flux to be deblended during the \ 2 calculation using 
the solved hydrogen underlying absorption, reddening, and temperature. In this way, the 
deblending itself affects the \ 2 an d therefore the relevant parameter solutions. The He (^ 889 ) 
flux is calculated by correcting the blended input for underlying absorption in H8, reddening 
the H8 intensity fraction, and subtracting it, 

^(3889) _ F(3889 + ff8)^(3889) + q g (ff8) . 046 mmc(m 

F(Hf3) F(Hf3) 1^(3889) 4 1 j 

(the same fraction is also subtracted from the equivalent width) . The error on this flux is not 
decreased so as to take the most conservative approach. Also, the errors on the equivalent 
widths, not being reported, are taken to be the larger of 10% or twice the fractional error in 
the flux. 



4.1. Investigating the Expanded Parameter Space 

Though the theoretical grounds for integrating the helium and hydrogen x 2 minimiza- 
tions are sound, this guarantees no practical benefit, and, in fact, the expansion of Monte 
Carlo and parameter space warns of more freedom for non-physical realizations. Synthetic 
testing provides some insight since the actual physical parameters are known. Characteristic 
values are chosen and used with the model equations to generate corresponding fluxes for 
all nine line ratios. For this analysis, n e = 100 cm -3 , a# e = 1.0 A, r = 1.0, T = 18,000 K, 
C(H/3) = 0.1, & H = 1-0 A, and £ = 1.0 x 10~ 4 , were used in conjunction with EW(H/3) = 250 
A and y+ = 0.08. Flux errors of 2% were assumed. The generating and solved parameters 
and abundances are listed in Table [2j 

Since the temperature and reddening are the only two parameters directly impacting 
both the helium and hydrogen x 2 , it might be expected that they would exhibit marked 
differences upon moving from sequential to integrated analysis. Each does as do their most 
strongly coupled parameters, the density and underlying hydrogen absorption. Integration 
improves the accuracy of the temperature marginally, but with a significantly larger un- 
certainty due to the added temperature dependence of the hydrogen lines. The reddening 
parameter, C(H/3) is essentially unchanged. The slight improvement for the reddening re- 
sults from its weaker sensitivity to perturbation due to its equal impact upon the helium 
and hydrogen lines. Since both the reddening and the hydrogen absorption affect the weaker 
hydrogen lines more strongly, they are tightly coupled such that a more accurate and precise 
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reddening correspondingly benefits the hydrogen absorption. The less accurate density with 
accompanying larger uncertainty stems from the previously discussed temperature-density 
degeneracy. The more uncertain temperature leads to the increased density uncertainty; 
while the non-linear nature of their trade-off biases the returned density. The benefit of 
the more accurate parameters propagates to the abundance determination, improving from 
0.08024 ± 0.00268 to 0.07995 ± 0.00324 (99.7% to 99.9%). The greater variance in the 
density and temperature, however, leads to an increased error bar. Overall, the effect of 
combining the helium and hydrogen calculations, on well-behaved spectra, is weak (and 
consistent at a set of higher and lower model parameters, though the sequential algorithm 
sometimes gives the marginally better abundance). The most notable feature is the similar- 
ity in the returned results, excepting a more pronounced effect of degeneracy in the density 
and temperature uncertainty. 

In the above synthetic analysis, both the sequential and integrated approaches achieve 
a well minimized hydrogen and helium x 2 ■ However, in real data, because the reddening 
and underlying hydrogen absorption both increase the flux ratio for the wavelengths shorter 
than H/3 (and both decrease |f§), minimization on the hydrogen lines alone may produce a 
relatively large x 2 ■ Then upon integration, the combination of the temperature dependent 
hydrogen emissivities, which were previously evaluated at the fixed [O III] temperature, and 
the density-temperature degeneracy may lead the derived temperature to deviate from the 
sequential value. The hydrogen x 2 exerts new pressure upon the temperature, but the weak 
(logarithmic) temperature dependence requires a large temperature adjustment. As a result, 
the temperature can slide substantially to reduce the hydrogen x 2 while maintaining the 
helium x 2 by inversely adjusting the density. The possibility for exactly this behavior is 



Table 2. Sequential and Integrated Synthetic Analysis 





Input 


Sequential 


Integrated 


Ho+/H+ 


0.08 


0.08024 ± 0.00268 


0.07995 ± 0.00324 




18,000 


18,212 ± 1768 


18,053 ± 3039 


N e 


100.0 


122.5 ± 123.4 


166.4 ± 307.9 


ABS(He I) 


1.0 


1.02 ± 0.08 


1.01 ± 0.10 


T 


1.0 


0.93 ± 0.40 


0.96 ± 0.60 


C(H/3) 


0.1 


0.095 ± 0.029 


0.096 ± 0.027 


ABS(H I) 


1.0 


1.34 ± 1.40 


1.19 ± 1.07 



a primary example of the necessity of accounting for systematic error via Monte Carlo in 
calculating the abundance error. 



4.2. Integrating the Dataset 

The resulting changes in the helium abundance in our nine objects due to the integration 
of the minimization procedure is shown in Figure [7J As evidenced in the preceding subsection, 
integration of the hydrogen and helium calculations could potentially yield significant changes 
to the temperature or reddening; however, for all ob jects, except one, the solved reddening 
was simil ar after in tegration to the result of IPaperTJ . The reddening for NGC 346 decreased 
from the iPaper it value of 0.174 to 0.156, which is a difference of about two sigma. This 



raises the helium abundance, as two of the three strong He lines are redder than H/3. 

The temperature changed significantly post integration for the three most metal poor 
objects (see Figure UJ). I Zw 18 showed a temperature increase (10,700 to 17,300 K) and 
a resulting increased abundance. This difference in temperature is significant, and results 
in an increase in the helium abundance of about 1.5 sigma. Note, however, that I Zw 18 
was not included in the original "high quality" dataset of Paper II because it did not meet 
two of the criteria (EW(H/3) and redshift), so this is likely due more to the quality of the 
spectrum and the looser constraint on the derived temperature. Nonetheless, the newly 
derived temperature is much more in line with that expected from the temperature derived 
from observations of the [O III] emission lines. 

However, the derived temperatures for SBS 0335-052 and SBS 940+544N both de- 
creased significantly (16,200 to 10,700 K and 18,500 to 12,800 K respectively) and showed 
corresponding helium abundance decreases. The newly derived temperatures are signifi- 
cantly different from those indicated by the [O III] emission lines. The magnitude of the 
effect on the abundance due to the temperature change for these three objects is apparent 
in Figure As discussed above, these shifting temperatures result from the combination of 
the temperature dependent hydrogen emissivities with the density-temperature degeneracy. 
Temperature changes allow the hydrogen \ 2 another degree of freedom for minimization, and 
the density compensates within the helium x 2 - This is supported by the fact that each of the 
above objects shows a larger contribution to the total \ 2 from the hydrogen lines than from 
the helium lines. Apart from inconsistency between the temperature and density used for the 
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hydrogen based corrections, integration would necessarily decrease the total x 2 ; i n validation 
of the procedure. The total x 2 of eight of nine objects does decrease (SBS 1159+545 increases 
from 1.35 to 1.70), but the large parameter variation witnessed is of concern. This provides 
further evidence for the importance of addressing the temperature-density degeneracy. 



5. Underlying Absorption and Wavelength Dependence 



Stellar absorption underlying the helium emission lines ha s long been recognized as po - 
tentially leading to und er estimates of th ei r mea sured fluxes. iDinerstein fe Shields! ( 119861 ); 



Izotov k Thuanl (I1998ah : ISkillman et all ((19981) and lOlive k Skillmanl (1200 if ) all demon 
strated examples of how failure to correct for underlying stellar absorption in the He I lines 
resulted in underestimates of nebular helium abundances. However, cor recting for the ef- 
fects of underlying He I a bsor ption is not straightforwar d (see discussions in 



2001 



Izotov et al 



Olive & Skillman 



20071 ). In lOlive k Skillmanl ( I200ll ). the simplest possible approach was 



suggested: namely, the assumption that the equivalent widths of all of the He I absorption 
lines were identical and could be obtained through the y 2 minimization for that value. The 
assumption of identical equivalent widths has merit since these lines are usually saturated 
in the stars that are producing them, and so the main problem is to essentially find a scal- 
ing factor to determine what fraction of the continuum is produced by the stars with these 
saturated lines. An alternative method is to use model stellar photospheres to produce an 
estimate of the relative strengths of the underlying absorption. In the second method, one 
needs estimates of the ages and metallicities of the stars and, in addition, the fraction of the 
continuum due to the stars. Note that the large equivalent widths of the Balmer lines in the 
"high quality" sample indicate that a significant fraction of the continuum is nebular, so the 
assumption of a purely stellar continuum is not entirely correct. 

There is evidence for two significant departures from the simplest assumption of iden- 
tical equivalent widths of underlying He I absorption. First, observations of stars indicate 
variations of order a factor of two (see discussion in 



Olive k Skillman 



20011 1 between dif- 



ferent He I lines, and this is also seen in the model stellar atmospheres. Second, different 
populations of stars contribute differently to different parts of the spectrum (i.e., young blue 
stars account for more of the blue continuum and older red stars account for more of the red 
continuum). If the stellar populations in the galaxies observed were single epoch bursts, then 
one could neglect the contributions to the continuum from different aged stellar populations. 
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Fig. 7. — Abundance comparison for the target objects in the sequential and integrated 
analyses. The red thinner lines are for the sequential analysis (using PFM emissivities and 
shown in Figure E]) with the integrated values given in black thicker bars. Again I Zw 18 
and NGC 346 are not used in calculating Y p . 
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However, it is becoming clear that the bursts of star formation that giv e rise to the bright H I 



regions typically studied have durations of order 100 Myr or longer (jMcQuinn et al.l 



2009 



and references therein). Since the He I absorption lines are formed in the bluer stars, the red 
continuum will have larger contributions from older populations, and thus, the equivalent 
widths of the redder He I lines will be weaker than predicted from models. 

In order to investigate the possible limitations of assuming identical equivalent widths 
of underlying stellar He I absorption, we constructed models of the emission from simple 

The code used 



stellar populations usin g the Violent Star Form a tion 



the input models from iGonzalez Delgado et al.l ( 120051 ) and 



egac y tool "SED(< 



>j~i 



Martins et al 



((2005J), and He I 



absorption line strengths as a function of stellar population age were examined. We then 
derived relative equivalent widths which were suitable over a large range of ages and these 
are presented here. We also did this for the H I absorption lines, although the variation is 
much smaller. 

The wavelength dependence of underlying absorption in terms of equivalent width for 
the hydrogen lines is normalized to H/3 while He A4471 is used for the helium lines. Those 
normalization values then represent the varying physical model parameters, a# and a# e , 

a{Ha) = 0.942 a(H/3) = 1.000 a{Hj) = 0.959 a(H5) = 0.896 
a(3889) = 1.400 a(4026) = 1.347 a(4471) = 1.000 
a(5876) = 0.874 a(6678) = 0.525 a(7065) = 0.400. (11) 

The value at H9 is also needed to successfully deblend He(A3889). A linear fit to the four 
hydrogen values above provides a value of 0.882. 

Our main goal here is to determine the relative importance of a secondary correction to 
the assumption of identical equivalent widths for underlying absorption. Thus, the precise 
choice of the wavelength dependence of the underlying absorption is not nearly as important 
as an assessment of the effect. 

Figure M shows the impact of assuming the wavelength dependent absorption given 
above. When comparing individual objects, in most cases there is no significant effect on the 
abundance. For most of the objects there is a very slight downward shift, presumably due to 
smaller corrections to the A5876 and A6678 lines. Since this effect is only a small correction 



'SED@ can be found at: |http://www. iaa.es/~mcs/secM 
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to the underlying absorption which itself is relatively weak, the minor impact is in agreement 
with expectations. Note that the impact is minor in large part due to the selection of targets 
for the "high quality" sample. By requiring large emission line equivalent widths for the 
targets, we ensure that this intrinsically uncertain correction will be small. 

Note that there are two objects where the difference caused by the new assumption 
is noticeable: SBS 0335-052 and NGC 346. In the case of SBS 0335-052, we note that the 
solution changes with almost every small change in assumptions that we test. This spectrum 
indicates a significant amount of self absorption in the helium lines, and, for large values of 
self-absorption our model is very uncertain. For NGC 346, the uncertainties are very small, 
so even small changes in assumptions have noticeable effects on the solution. 



6. Neutral Hydrogen Collisional Emission 

To correct the measured hydrogen fluxes for emission due to neutral hydrogen collisional 
excitation, the relative amount of collisional to recombination emission for a given line can 
be calculated as follows, 

C {x) = n(H I^K^BR^BR^ = n(H I) K eff = 

R n(H II)a + ^jBRj^2n e n(H II) a e /f a e ff ' 

where K represents the collisional transition rate from the ground state to some level, i, 
above the transition level of interest, j (j — > 2 = A). Downward transitions from i occur 
with a variety of branching paths so BR^j represents the relevant branching fraction for 
transitioning ultimately to j. a is the effective recombination rate to level j {ot+^j = u e ff 
with "+ — > j" serving as an analogue to i^i-n but for ionization recombining to level j). 
K e ff is the sum over the excitation levels of the collisional transition rate, i^i-^, times the 
corresponding branching fraction, BR^j. Given the relative complexity of the definition 
of ^r, it is useful to note that the electron density and branching fraction from level j to 2 
cancel, and we are left with the ratio of the neutral hydrogen density, n(H I), to the ionized 
hydrogen density, n(H II), times a collisional rate over a recombination rate. The ratio 
n(H ii) ^ S then defined as the model parameter £ which will be solved via the minimization. 

As defined above, ^ makes the simplifying but very accurate assumption, at these tem- 
peratures, that all of the neutral hydrogen is excited from the ground state. The collisional 
rates are based on power law fits (AT B ) to the effective collisional strengths, T, reported 
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Fig. 8. — Abundance comparison for the target objects upon inclusion of wavelength depen- 
dent absorption. The red thinner lines are for the integrated analysis (of Figure [7]) with the 
same analysis performed including wavelength dependent absorption given in black thicker 
bars. As before, I Zw 18 and NGC 346 are not used in calculating Y p . 
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ni 



Anderson et al 



( 120021 ) over the temperature range 5000 to 35,000 K. The effective colli- 
sional strength as defined below allows for more convenient fitting due to its more gradual 
temperature dependence and is easily converted to the collisional excitation rate via the 
accompanying rate formula. 

f°° e ■ 

T« = / Q(i i) expf-^)d(- 



13 



k B T knT 



13 ^ 



(13) 



Ui V k B T k b T 



/ 1 __,-13.6(l-ir) 



Km = 4.004 x — «p( hT ^ )T M 

with fi, the collisional strength, related to the cross-section by Qjj = 2me ^ i aij(e i ), ej is the 
final outgoing electron energy, e, is the incident electron energy, In = —13.6 eV the ionization 
potential, a is the Bohr radius, and = (2s + 1)(2Z + 1) the statistical weight of the level. 
The branching fractions are calculated directly from the Einstein transition coefficients, and 
the effective recombination rate is fit from the case B data of Hummer Sz Story (1987) for an 
electron density of 100 per cm 3 and over the temperature range 10,000 to 30,000 K. Again a 
power law fit is chosen for simplicity and models well. In theory, the collisional sum includes 
an infinite number of levels, but the probabilities fall off rather quickly. For this work the 
sum excludes any terms contributing less than 1%, giving the following contributions: 
For Ha, 3s, 3p, 3d, 4s, 4d, and 4f are included; 
For H/3, 4s, 4p, 4d, 4f, 5s, 5d, 5f, and 5g; 
For H7, 5s, 5p, 5d , 5f, and 5g; 



Anderson et al 



(120021 ) only report values for collisions up to a principle quantum number of 
5; therefore, to calculate the contribution to H<5, the equation for H7 was used and scaled 
by the appropriate energy level difference, 

R {6) = R il)eM k-f } - (15) 

Table [3] lists the coefficients for constructing ^ £i£i , the quantity that multiplies £ and 

T \ 



provides the increasing collisional enhancement factor for increasing temperature (T 4 — j^j- 

;>>exp(-^)T 4 c « (16) 



K eff V- , bi 

at ef f ^— ' 1 A 
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As can be clearly seen in Figure M below, the collisional excitation grows strongly with 
temperature. This has the effect of leaving £ poorly constrained at low temperature, but 
potentially well constrained at high temperature due to its now significant effect on the 
hydrogen x 2 - The neutral hydrogen collisional correction has no direct effect on the helium 
X 2 ', however, it will directly raise the abundance due to the decrease in the normalizing 
H/5 flux. In this regard, when working with lower temperature objects (like NGC 346) the 
program is vulnerable to producing non-physical results (very high neutral fractions). Here 
we will analyze NGC 346 for demonstration purposes, but we do not use NGC 346 in our 
final determination of the primordial helium abundance. 



6.1. A Return to Integration 

The effects of integrating the helium and hydrogen \ 2 minimizations were already dis- 
cussed in section 14.11 However, apart from the general motivation of increased consistency, 
the inclusion of the neutral hydrogen fraction provides a strong practical benefit. Again, 
synthetic testing, with its known parameters, provides insight. The same generating param- 
eters as in section 1431 are used with the addition of £ = 1.0 x 10~ 4 . The generating and 
solved parameters and abundance are listed in Table HI 

As before, the most significant changes are seen in the temperature, density, reddening, 
and underlying hydrogen absorption. Now, the addition of neutral hydrogen requires solving 
three parameters (C(H/3), a#, £) from the hydrogen lines alone if sequential analysis is used, 
but only two for integrated analysis (since the reddening affects both the hydrogen and helium 
lines). This dramatically improves the reddening accuracy, 52% to 82%. The aforementioned 
similarity of wavelength behavior between the reddening and hydrogen absorption again 
leads to corresponding improvement for the hydrogen absorption. The much larger and 
more discrepant neutral hydrogen fraction, with an especially increased uncertainty, is, on 
the surface, disconcerting. The increased temperature uncertainty is the cause. As discussed 
in the previous section, at low temperatures the neutral hydrogen collisional emission rate is 
greatly reduced permitting a larger fraction of neutral hydrogen. Therefore the introduction 
of lower temperature solutions after integration leads to lower neutral hydrogen collisional 
emission rate which is then compensated by a higher fraction of neutral hydrogen. The 
resulting correction is actually smaller for the integrated analysis (e.g., £j^(H/3) = 3.3% — > 
1.4%) such that the neutral hydrogen correction is improved from the sequential to integrated 
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Table 3. Coefficients for the neutral hydrogen collisional correction 



Line: 


Coefficient 


Terms 
















Ha 


a 
b 

c 


0.4155 
-14.80 
0.4209 


2.4965 
-14.03 
0.5853 


2.4063 
-14.03 
0.6187 


0.2914 
-14.80 
0.6766 


0.3685 
-14.80 
0.7076 


4.6426 
-14.03 
0.7788 






H/3 


a 
b 

c 


0.2384 
-15.15 
0.3082 


0.6964 
-14.80 
0.4978 


0.1991 
-15.15 
0.6017 


0.1409 
-15.15 
0.6765 


0.2201 
-15.15 
0.7293 


1.9228 
-14.80 
0.7535 


1.4845 
-14.80 
0.7845 


2.8179 
-14.80 
0.9352 


H7 


a 
b 

c 


0.3629 
-15.15 
0.3598 


0.8351 
-15.15 
0.6533 


2.0044 
-15.15 
0.7281 


1.4757 
-15.15 
0.7809 


2.7947 
-15.15 
0.8582 








m 


a 
b 

c 


0.3629 
-15.34 
0.3598 


0.8351 
-15.34 
0.6533 


2.0044 
-15.34 
0.7281 


1.4757 
-15.34 
0.7809 


2.7947 
-15.34 
0.8582 









Table 4. Comparing Sequential and Integrated Analyses with Synthetic Data 





Input 


Sequential 


Integrated 


Hc+/H+ 


0.08 


0.08375 ± 0.00502 


0.08099 ± 0.00373 


T e 


18,000 


17,138 ± 1812 


17,505 ± 3103 


N e 


100.0 


141.9 ± 252.7 


173.7 ± 311.4 


ABS(He I) 


1.0 


1.16 ± 0.10 


1.05 ± 0.10 


T 


1.0 


1.35 ± 0.42 


1.15 ± 0.68 


C(H/3) 


0.1 


0.05 ± 0.05 


0.08 ± 0.03 


ABS(H I) 


1.0 


2.16 ± 1.65 


1.32 ± 1.10 


£ x 10 4 


1.0 


11.62 ± 12.22 


60.65 ± 269.67 
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Fig. 9. — The collisional correction factor for four hydrogen emission lines due to neutral 
hydrogen. The behavior is dominantly exponential with temperature, with, as expected, 
larger corrections for the lower energy emission lines. Note that the enhancement will scale 
linearly with increasing the neutral to ionized hydrogen ratio, £ (a characteristic value of 
10~ 4 is used for the plot). 
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analysis. The full benefit of the more accurate solution is manifested in the abundance, 
improving from 0.08375 ± 0.00502 to 0.08099 ± 0.00373 (95% to 99%). 

To summarize, the use of the solved temperature in evaluating the hydrogen lines, 
though necessary for consistency, exposes the neutral hydrogen fraction to the temperature- 
density degeneracy, as evidenced by the uncertainty on the neutral hydrogen fraction. Apart 
from this temperature feedback, however, integration clearly benefits the neutral hydrogen 
solution as well as the other parameters. Primarily, the reddening, constrained through the 
helium and hydrogen lines, is more precisely determined; therefore, the underlying absorption 
and neutral hydrogen fraction are allowed less range in correcting the hydrogen lines, and 
the abundance determination benefits. 

6.2. The Dataset Analysis with the New Neutral Hydrogen Parameter 

The introduction of the neutral hydrogen collisional emission complicates the abundance 
determination. A three way degeneracy between electron density, temperature, and the 
neutral hydrogen fraction emerges. These degeneracies are shown in Figure [10] which shows 
the Monte Carlo results for NGC 346. Plotted are the individual results for n e ,T, and £ 
over the 1000 MC realizations. As one can see, the solutions map out a wide range for these 
parameter values. 

Since the neutral hydrogen fraction multiplies — directly, the neutral hydrogen fraction 
can be increased to compensate for a decrease in the collisional to recombination rate which 
is itself only temperature dependent. Furthermore the neutral hydrogen collisional emission 
correction can have a surprisingly significant effect in raising the abundance as evidenced 
in Figure [TTJ Since the collisional correction for H/3 multiplies each helium line abundance, 
it multiplicatively raises the abundance but does not change the value of the corresponding 
X 2 - As a result, the abundance of all nine objects increases as does the derived primordial 
helium mass fraction. 

6.3. Comparing Implementations 

The approach detailed above differs from the recent analogous additions to ITS07 and 
PLP07 in several important ways. ITS07 varies the collisional contribution to Ha between 
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Fig. 10. — Monte Carlo plots (1000 points) of density, temperature, and neutral hydrogen 
fraction for NGC 346. The strong correlation between the three parameters demonstrates 
the difficulty in solving for them independently. 
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Fig. 11. — Abundance comparison for the target objects detailing the effect of the neutral 
hydrogen collisional emission. The red thinner lines are for the integrated analysis including 
wavelength dependent absorption (of Figure [H]) with the same analysis performed including 
neutral hydrogen collisional emission given in black thicker bars. As before, I Zw 18 and 
NGC 346 are not used in calculating Y p . 
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and 5% and takes the collisional contribution to H/3, and therefore to y + , to be | as large. 
The minimum of the reddening corrected helium \ 2 (each minimized over n e , r, and T) is 
taken to be the solution. This method does not allow for a full sampling of possible parameter 
values as we find in the MC analysis exhibited in Figure [10J Regarding the neutral hydrogen 
collisional emission, the key differences with this work are modeling the relative collisional 
contributions to Ha and H/3 only in determining the reddening and treating the collisional 
emission fraction in its entirety as a model parameter rather than a neutral hydrogen fraction 
times a temperature dependent relative collisional emission rate. 

PLP07 creates a photoionization model for each of their objects using CLOUDY which 
outputs the collisional contributions to the hydrogen lines. Their model also makes use of 
the collisional strength data of Anderson et al. (2002). The reddening is calculated from 
the corrected Ha and H/3, and y + is increased by the decrease in H/3. The primary dif- 
ferences with the neutral hydrogen emission of this work are, again, correcting only Ha 
and H/3 in determining the reddening and the construction of a photoionization model, in 
particular requiring the specification of ionizing radiation and geometry of the region, to de- 
termine the neutral hydrogen collisional contribution. These differences extend beyond the 
neutral hydrogen collisional contribution modeling and determination. We emphasize that, 
here, in addition to the collisional contributions to H7 and US, all six He lines, due to the 
combined hydrogen and helium calculations, are included in determining the reddening as 
well. Furthermore, the integrated approach utilizes the same temperature in determining the 
hydrogen recombination emission, hydrogen collisional emission, and the helium emission. 

Table [5] below details the relative collisional to recombination emission contributions, 
^, for the objects here analyzed with corresponding corrections of ITS07 and PLP07 listed 
where available. Note that for this work ^ is computed for each Monte Carlo realization 
from the solved neutral hydrogen fraction, £, and temperature, T, with the table listing the 
Monte Carlo average. The magnitude of the corrections calculated for this work is in broad 
agreement with that found in ITS07 (for which the solved Ha correction is restricted to less 
than 0.05) and PLP07. The results for individual objects, however, do not in general agree. 
The most notable disagreements are for Haro 29 and NGC 2363A, where this work finds a 
much smaller correction than either ITS07 or PLP07, and for Mrk 193, much larger than 
ITS07. Close agreement is found for Ha for SBS 0335-052 with PLP07 and SBS 0940+544N, 
SBS 1159+545, and SBS 1420+544 (only reported by ITS07). As evidenced by those last 
three objects, the temperature dependent neutral hydrogen collisional contribution of H/3 
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relative to Ha is always less than |, the value used by ITS07. PLP07 finds this ratio of 
the H/3 to Ha contribution always greater than this work but also always less than ITS07. 
Both this work and PLP07 calculate relative corrections to the Ha emission in excess of 
the 5% limit imposed by ITS07 for several objects. Overall, the different neutral hydrogen 
emission models yield qualitatively and quantitatively different results, but the variance and 
magnitude of the effect are in accordance. 



Table 5. Comparison of Neutral Hydrogen Collisional Contribution 





Ha 




H/3 




H7 




H<5 




Ha 


H/3 


Ha 


H/3 


Object 












This Analysis 












ITS07 


PLP07 a 


I Zw 18 


0.0070 


± 


0.0089 


0.0037 


± 


0.0083 


0.0030 


± 


0.0083 


0.0027 


± 


0.0083 


0.0186 


0.0222 


0.075 


0.056 


SBS 0335-052 


0.0881 


± 


0.0607 


0.0425 


± 


0.0464 


0.0334 


± 


0.0488 


0.0297 


± 


0.0501 


0.0033 


0.0036 


0.094 


0.071 


SBS 0940+544N 


0.0545 


± 


0.0323 


0.0273 


± 


0.0223 


0.0217 


± 


0.0237 


0.0194 


± 


0.0247 


0.0499 


0.0490 






SBS 1159+545 


0.0542 


± 


0.0342 


0.0262 


± 


0.0244 


0.0205 


± 


0.0259 


0.0182 


± 


0.0268 


0.0492 


0.0452 






SBS 1420+544 


0.0579 


± 


0.0373 


0.0309 


± 


0.0281 


0.0249 


± 


0.0293 


0.0226 


± 


0.0301 


0.0497 


0.0520 






Haro 29 


0.0011 


± 


0.0016 


0.0006 


± 


0.0015 


0.0004 


± 


0.0016 


0.0004 


± 


0.0016 


0.0357 


0.0426 


0.034 


0.021 


Mrk 193 


0.0442 


± 


0.0400 


0.0191 


± 


0.0339 


0.0145 


± 


0.0348 


0.0126 


± 


0.0353 


0.0013 


0.0004 






NGC 2363A 


0.0013 


± 


0.0069 


0.0005 


± 


0.0069 


0.0004 


± 


0.0069 


0.0004 


± 


0.0069 


0.0115 


0.0156 


0.038 


0.028 


NGC 346 


0.0461 


± 


0.0282 


0.0193 


± 


0.0164 


0.0145 


± 


0.0181 


0.0125 


± 


0.0191 






0.011 


0.007 



a In terms of the variable used in PLP07, x x = {^l"" 1 , §(A) = - 

I(\)tot S> ' 1 
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Ionization Correction Factors 



If the ionized hydrogen is not spatially co-existent with the ionized helium, then a 
correction needs to be made for this effect. This can happen for two different reasons. If the 
ionizing radiation field lacks sufficient high energy photons, then some of the helium will be 
neutral where the hydrogen is ionized. It is also possible that lower energy ionizing photons 
will be preferentially captured by hydrogen atoms resulting in regions where helium is ionized 
and hydrogen is not (this can occur at the edges of 1-dimensional photoionization models 
of H II regions). Since neutral inclusions are impossible to detect in spectra of unresolved 
H II regions, corrections for these effects are inherently uncertain. Typically photoionization 
models are used as a guideline for the importance of this correction. 

The consensus approach to dealing with neutral helium has been to observe objects 
with spectra that indicate that the pre s ence of neutral helium is highly unlikely. Based 
on photoionization models, iPagel et al.l ( 119921 ) established a criterion to select only H II 
regions with sufficiently hard radiation fields, and this method has been followed through 
most investigations. 

Th e reverse correction, dealing with the pre s ence of neutral hydrog en, has been dis- 



cussed (iViegas et al. 



2000 



Ballantyne et al 



2000 



Gruenwald et al 



200 2]), but gene rally is 



not applied nor accounted for in the uncertainty. Interestingly, iBallantyne et al.l (120001 ) found 
that by eliminating H II regions with relatively strong [O I] A6300 emission, that the scatter 
in helium abundances decreased, and interpreted this as a sign of neutral hydrogen within 
the ionized helium. It is also possible that the [O I] A6300 emi ssion is indicati ve of collisional 



excitation due to the presence of supernova remnants (e.g., ISkillmanl Il985l ). and that the 
elevated A6300 emission could simply be an indicator that the nebula is not excited purely 
by photoionization, and thus the simple model used to derive abundances is not appropriate. 

In this regard, detections of collisionally excited Balmer emission are of interest. If as 
much as 1% of the helium is co-existent with neutral hydrogen, then collisional excitation of 
the lower Balmer emission lines would be important in that volume. To our knowledge, none 
of the studies of these ionization correction factors ever looked into detecting the neutral 
hydrogen by looking at the collisionally excited Balmer emission. Although it is beyond the 
scope of this work, we believe that this would be a worthwhile pursuit. For the present, 
we estimate that the amplitude and the uncertainty in the ionization correction are much 
smaller than the other factors discussed in this paper. 
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8. Primordial Helium 



Figures El El EH summarize the results of each s ubsequent change or addition of 
sections [2H6J Table [6] provides numerical results comparing Paper I J and the analysis of this 
work with all of the changes and addition s included. The abundance and all seven (except 
th e neutral hydrogen fraction for lPaper III ) parameters are listed. I Zw 18 was not analyzed 



in 



Paper III and is therefore not listed. The cumulative effect of the updated emissivities, 



integration, wavelength dependent absorption, and neutral hydrogen collisional emission is to 
raise the abundance for seven of the eight objects. SBS 1159+545 undergoes the previously 
mentioned decrease from the updated emissivities but is then very stable. The abundance 
error increases for six of the eight with Haro 29 and NGC 2363A stable after integration 
decreased their error. Mutual objects with ITS07 and PLP07 allows for further comparison 
(see Table Cj). Broadly, this work finds lower temperatures, higher densities, and lower optical 
depths than ITS07. There is no clear pattern with the abundance though on average it is 
higher. Our derived densities are lower than those of PLP07 with correspondingly higher 
abundances. The only object not in agreement is NGC 346. The solution for NGC 346 
exhibits susceptibility to large variance, relative to its small uncertainty. 

Our primary interest is the primordial helium abundance (mass fraction), Y p . To ex- 
tr apolate ba ck to a primordial value, the oxygen to hydrogen mass fraction, O/H, calculated 



in 



Paper III , is used in conjunction with Y, the heli um mass fraction, including a correction 
for y^E The relevant values are listed in Table El Paper U analyzed NGC 346 but did not 



include it in the regression because its more evolved chemical nature invokes mor e error in 
assuming a linear relationship between y + and O/H. I Zw 18 was not analyzed in iPaper II 
both because of its low equivalent width for H/3 (EW(H(3) = 135 < 200A) making it suscep- 
tible to large corrections for underlying absorption and because of it's radial velocity making 
it susceptible to contamin ation by Galactic Na I absorption. A linear regression of Y for the 
seven targets of lPaper II] versus O/H yields, 



Y p = 0.2561 ±0.0108, 

with slope 9 ± 192 and x 2 = 2.0. This result is shown in Figure 
uncertainty, this agrees well with the WMAP result of Y p = 0.2486 ± 0.0002. IPaper II 



(17) 

Given the large 



found 



5 This work takes Z = 20(O/H) such that Y = MlzMg/gl) 
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Fig. 12. — Helium abundance (mass fraction) versus oxygen to hydrogen ratio regression 
calculating the primordial helium abundance. 
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Y p = 0.2495 ± 0.0092. The increased intercept and error of this work is directly attributable 
to y + increasi ng for six targets and an uncertainty increase for four. The slope is similar to 
that found in lPaper III (54 ± 187) and still consistent with zero. Given the short baseline, a 



mean evaluation yields, 



Y p = 0.2566 ± 0.0028 



for seven which is very similar to the Bayesian result (IHogan et al.lll997l ) of 0.2563 ± 0.0028. 



Because of the large uncertainty in Y, Equation f)17p is in agreement with PLP07, 
Y p = 0.2477 ± 0.0029, using PFM, and ITS07, Y p = 0.2516 ± 0.0011, using BSS rescaled to 
PFM emissivities. The striking, larger error found here (as in Paper II) in comparison to 
PLP07 and ITS07 is attributable to use of an unrestricted Monte Carlo, the application of 
all systematic corrections simultaneously, and a smaller population of objects (than ITS07). 
I Zw 18 is below the regression line and, being the lowest metallicity point, it would lower 
the intercept and raise the slope. NGC 346 would extended the O/H domain and decrease 
the uncertainty in both the intercept and slope. If both regions are included, we would find, 
Y p = 0.2528 ± 0.0060, with a slope of 64 ± 80. 

As noted earlier, the inclusion of the neutral hydrogen collisional correction has a large 
impact on the result. As seen in Figure [TT| the same seven object regression performed using 
the analysis with £ = yields a lower tighter intercept of, Y p = 0.2394 ± 0.0083. Though 
not desirable, the spread of the results with and without the neutral hydrogen correction 
highlights the dominance of systematic effects and therefore the necessity of quantifying the 
full effect the associated error. 

Since the main goal of this paper is an improved assessment of the errors on the helium 
abundances derived from spectra of metal poor H II regions, it is useful to compare the sizes of 
the error bars determined in this study with those determined from the same spectra in other 
studies. Figure [13] shows the sizes of the total errors associated with each helium abundance 
determination for all of the objects studied here. One surprising result of this comparison is 
that the relative size of the errors is a strong function of oxygen abundance for all studies. 
To our knowledge, this has never been seen or discussed before. One possible explanation 
is that the most metal poor objects with suitably bright H II regions are relatively rare, so 
that the quality of the spectra for these objects may be lower. Note also that typically the 
errors that we derive for each spectrum are larger than previously derived by others. The 
errors derived by ITS07 are typically the smallest, and this is directly attributable to their 
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13. 



A comparison of the c alculated uncerta i nties on the ind i vidual 



data points from three different st udies ([Olive fc Skillmanl 12004 llzotov et all 120071 ; 
Peimbert. Luridiana. fc Peimbertl 120071 and this work). The data points are ordered from 
lowest oxygen abundance to highest oxygen abundance. Note the trend of decreasing uncer- 
tainty with increasing abundance. Typically our uncertainties are larger than other previous 
analyses of the same spectra. 
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method of deriving the electron temperature which limits that value to a relatively small 
range. PLP07 derive electron temperatures in a similar method to ours, but do not include 
a Monte Carlo analysis of the errors. Thus, the errors of PLP07 are typically intermediate 
to those of ITS07 and those presented here. 



Table 6. Comparison of Physical Conditions and He + /H + Abundance Solutions 



Object 


He+/H+ 




N e 


ABS(He I) 


T 


C(H/3) 


ABS(H I) 


£ x 10 4 




^arjei^JI 


SBS 0335-052 


0.0763 ± 0.0049 


15,940 ± 2710 


347 ± 942 


0.1 ± 0.2 


4.6 ± 0.9 


0.121 ± 0.014 


1.3 ± 1.9 




SBS 0940+544N 


0.0841 ± 0.0035 


19,260 ± 2480 


35 ± 140 


0.1 ± 0.3 


0.4 ± 0.4 


0.048 ± 0.023 


0.0 ± 1.4 




SBS 1159+545 


0.0838 ± 0.0031 


19,330 ± 2100 


75 ± 116 


0.0 ± 0.1 


0.5 ± 0.4 


0.065 ± 0.016 


0.2 ± 1.1 




SBS 1420+544 


0.0854 ± 0.0026 


20,370 ± 1740 


27 ± 83 


0.12 ± 0.18 


3.32 ± 0.35 


0.16 ± 0.03 


0.4 ± 0.9 




Haro 29 


0.0844 ± 0.0024 


16,680 ± 890 


33 ± 97 


0.5 ± 0.2 


0.4 ± 0.2 


0.001 ± 0.004 


2.7 ± 0.3 




Mrk 193 


0.0845 ± 0.0023 


15,643 ± 1130 


23 ± 95 


0.37 ± 0.14 


2.29 ± 0.35 


0.254 ± 0.019 


0.36 ± 0.51 




NGC 2363A 


0.0801 ± 0.0046 


14,040 ± 1060 


285 ± 343 


0.35 ± 0.28 


1.4 ± 0.4 


0.117 ± 0.003 


1.4 ± 0.3 




NGC 346 


0.0828 ± 0.0008 


13,420 ± 280 


2 + 19 


0.09 ± 0.09 


0.01 ± 0.04 


0.174 ± 0.008 


0.0 ± 0.5 




Re- analysis 


I Zw 18 


0.08008 ± 0.00469 


17,588 ± 2951 


76.7 ± 277.2 


0.42 ± 0.36 


0.42 ± 0.56 


0.00 ± 0.00 


4.76 ± 0.82 


7.89 ± 70.34 


SBS 0335-052 


0.08536 ± 0.00644 


15,822 ± 2395 


120.8 ± 254.4 


0.39 ± 0.31 


5.32 ± 0.87 


0.05 ± 0.04 


1.82 ± 1.04 


117.44 ± 262.19 


SBS 0940+544N 


0.08506 ± 0.00484 


16,689 ± 2124 


149.3 ± 277.0 


0.09 ± 0.20 


0.61 ± 0.56 


0.00 ± 0.00 


0.03 ± 0.17 


25.82 ± 33.46 


SBS 1159+545 


0.08332 ± 0.00390 


16,177 ± 2318 


325.2 ± 330.0 


0.03 ± 0.06 


0.39 ± 0.43 


0.02 ± 0.02 


0.59 ± 0.63 


68.26 ± 152.42 


SBS 1420+544 


0.08939 ± 0.00393 


20,214 ± 2695 


83.1 ± 234.0 


0.18 ± 0.16 


2.89 ± 0.62 


0.16 ± 0.03 


0.02 ± 0.14 


30.30 ± 200.16 


Haro 29 


0.08576 ± 0.00174 


17,092 ± 665 


37.1 ± 48.3 


0.52 ± 0.15 


0.13 ± 0.16 


0.00 ± 0.00 


3.38 ± 0.24 


0.20 ± 0.39 


Mrk 193 


0.08660 ± 0.00435 


13,924 ± 2219 


117.1 ± 319.9 


0.45 ± 0.19 


2.63 ± 0.78 


0.21 ± 0.04 


0.64 ± 0.74 


1003.06 ± 8465.55 


NGC 2363A 


0.08328 ± 0.00244 


13,578 ± 992 


194.2 ± 203.1 


0.46 ± 0.18 


1.55 ± 0.34 


0.12 ± 0.01 


0.29 ± 0.44 


5.20 ± 49.25 


NGC 346 


0.08777 ± 0.00168 


12,847 ± 708 


30.2 ± 70.2 


0.37 ± 0.15 


0.03 ± 0.09 


0.13 ± 0.01 


0.10 ± 0.28 


358.72 ± 758.49 
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9. Discussion & Future Improvements 

The preceding work in Papers I and II emphasized the importance of the systematic 
uncertainties in calculating y + and the need for Monte Carlo analysis. Extending that 
work to allow for self-consistent analysis reflecting the common temperature and reddening 
parameters between the hydrogen and helium lines does not diminish the importance of 
the systematic uncertainties or the case for the larger errors determined via Monte Carlo. 
In fact, the addition of the neutral hydrogen collisional emission, introduced with a new 
physical parameter, the neutral hydrogen fraction, highlights that abundance errors may be 
even larger than previously thought. In addition to the need for observations to constrain 
systematic uncertainty as much as possible, sections [3] & 0] both demonstrate the significant 
impact that relatively small systematic errors in the observed fluxes, even of the weakest 
lines, can have on y + . Clearly high quality spectra are critical for meaningful analysis. 

Given higher resolution spectra, the emission line fluxes could be measured separately 
from the effects of underlying absorption. This is possible due to collisional broadening of 
the stellar absorption lines, such that the emission line can be measured relative to its local 
continuum. In other words, the absorption effect in the continuum would be visible in depth 
and extent with a narrow emission line at the minimum. As a result the flux could be mea- 
sured relative to the absorption feature rather than the broader continuum. The underlying 
helium and hydrogen (through H/3) lines both directly impact the abundance determination. 
Furthermore, decreasing the number of solution parameters (by eliminating the parameters 
associated with underlying absorption) will in general decrease the uncertainty. Therefore, 
an improvement in the abundance determination is ensured. Typically, underlying Balmer 
emission is much easier to detect in the spectrum of a metal poor H II region than underlying 
He I emission, so starting with eliminating the correction for underlying Balmer absorption 
is the most promising. 

Though difficult to gauge theoretically, some insight into the degree of improvement from 
more accurate measurements can be gained through synthetic data. Assuming a typical value 
of 1.0 A for the underlying Hydrogen and Helium absorption, we can generate corresponding 
fluxes when the remaining parameters are chosen. For r = 1, C(H/3) = 0.1, T = 18000 K, 
£ = 10 -4 , and y + = 0.08, we generated fluxes to determine the helium abundance when 
either the underlying absorption is assumed to be known or unknown. At a low density of 
10 cm -3 there was no improvement in the determination of y + due to a large bias in the 
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Table 7. Comparison with ITS07 & PLP07 



Analysis 



C(H/3) 



O/H 



I Zw 18 



ITS07 


0.2485 ± 0.0081 


18,100 ± 500 


11 +™ 

11 -0 


1 o +0.00 
l.OU _g gg 




0.16 ± 0.01 




0.2415 ± 0.0104 


18,600 ± 600 


132 t\f 


42 +0 ' 51 
U - 4Z -0.25 




0.15 ± 0.01 


PLP07 


0.2505 ± 0.0087 




90 ± 80 


0.06 ± 0.05 


0.02 ± 0.01 


0.22 ± 0.03 
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Monte Carlo stemming from the requirement that the density remain positive. At a more 
characteristic density near 100 cm -3 the benefit was significant. The returned value for the 
helium abundance was y+ = 0.080997 ± 0.003735 and y + = 0.080002 ± 0.002841 for the 
underlying absorption solved and removed cases respectively. 

Ultimately, the sensitivity of the determination of the primordial helium abundance de- 
mands that analysis be self-consistent and comprehensive; unfortunately efforts in providing 
such an analysis yield uncertainties that limit the utility of that determination. Nevertheless, 
fully understanding the terrain of potential obstacles is necessary to direct further measure- 
ment efforts so as to provide a reliable calculation of the primordial helium abundance. 
Apart from improving the spectra, the neutral hydrogen collisional correction for the weaker 
hydrogen lines can be easily improved from extended principal quantum number collisional 
strengths. Furthermore the use of a further hydrogen line would benefit the neutral hydro- 
gen fraction and hydrogen absorption since those two parameters are determined exclusively 
through the three hydrogen lines. The next Balmer line, He (A3970) is blended with Ne III 
(A3967). The contribution of Ne III A3967 to the He A3970 could be corrected for by scaling 
from the stronger Ne III A 3868 line, or by observing at sufficiently higher spectral resolution 
in order to deblend the two lines. H8 (A3889), is blended with He I A3889, and would not 
be helpful independent for this type of analysis. This leaves the relatively very weak H9 
(A3835) and H10 (A3798) as the best candidates. The use of other emission lines which are 
sensitive to the density may also be useful in reducing the uncertainty due to the degeneracy 
between density and temperature. However, most candidates with the appropriate electronic 
structure are not the dominant ionic species in metal-poor H II regions (e.g., [S II], [O II]). 
A possible exception is the [CI III] emission line pair AA5518, 5538. 

In summary, it should be emphasized that even though the precision of the primordial 
helium abundance is not increased, the benefit of the analysis improvements is precisely in 
demonstrating the full effect of the systematic errors in restricting the abundance determi- 
nation. 

Note Added : On completion of this paper, we became aware of the newest result of 



Izotov fc Thuan! (120101 ) . While we find a very similar central value for Y p , we disagree with 
their conclusion calling for physics beyond the standard model of BBN. The cumulative 
uncertainties described above are sufficiently large that, within these uncertainties, standard 
BBN is not presently threatened on the basis of 4 He observations. 



-47- 



We would like to Yuri Izotov for providing additional information concerning his spectra. 
Preliminary results of this work were presented at the conference IAU Symposium 268, "Light 
Elements in the Universe," and we would especially like to thank the conference organizers 
Corinne Charbonnel and Monica Tosi for the opportunity to do so. We would like to thank 
Y. Izotov, A. Peimbert, M. Peimbert, G. Steigman, E. Terlevich, and R. Terlevich, for 
informative and valuable discussions. The work of KAO is supported in part by DOE grant 
DE-FG02-94ER-40823. EDS is grateful for partial support from the University of Minnesota. 
This research has made use of NASA's Astrophysics Data System Bibliographic Services and 
the NASA/IPAC Extragalactic Database (NED), which is operated by the Jet Propulsion 
Laboratory, California Institute of Technology, under contract with the National Aeronautics 
and Space Administration. 



-48- 
REFERENCES 

Anderson, H., Ballance, C. P., Badnell, N. R., & Summers, H. P. 2002, Journal of Physics 
B Atomic Molecular Physics, 35, 1613 

Ballantyne, D. R., Ferland, G. J., & Martin, P. G. 2000, ApJ, 536, 773 

Benjamin, R. A., Skillman, E. D., & Smits, D. P., 1999, ApJ, 514, 307 

Benjamin, R. A., Skillman, E. D., & Smits, D. P. 2002, ApJ, 569, 288 

Cardelli, J. A., Clayton, G. C, & Mathis, J. S. 1989, ApJ, 345, 245 

Coc, A., Vangioni-Flam, E., Casse, M., k Rabiet, M. 2002, Phys. Rev. D, 65, 043510 

Coc, A., Vangioni-Flam, E., Descouvemont, P., Adahchour, A., & Angulo, C. 2004, ApJ, 
600, 544 

Cuoco, A., Iocco, F., Mangano, G., Miele, G., Pisanti, O., & Serpico, P. D. 2004, Interna- 
tional Journal of Modern Physics A, 19, 4431 

Cyburt, R. H. 2004, Phys. Rev. D, 70, 023505 

Cyburt, R. H., Fields, B. D., & Olive, K. A. 2001, New Astronomy, 6, 215 

Cyburt, R. H., Fields, B. D., & Olive, K. A. 2002, Astroparticle Physics, 17, 87 

Cyburt, R. H., Fields, B. D., & Olive, K. A. 2003, Physics Letters B, 567, 227 

Cyburt, R. H., Fields, B. D., & Olive, K. A. 2008, Journal of Cosmology and Astro-Particle 
Physics, 11, 12 

Cyburt, R. H., Fields, B. D., Olive, K. A., & Skillman, E. 2005, Astroparticle Physics, 23, 
313 

Davidson, K., & Kinman, T. D. 1985, ApJS, 58, 321 

Descouvemont, P., Adahchour, A., Angulo, C, Coc, A., & Vangioni-Flam, E. 2004, Atomic 
Data and Nuclear Data Tables, 88, 203 

Dinerstein, H. L., & Shields, G. A. 1986, ApJ, 311, 45 



Esteban, C, Bresolin, F., Peimbert, M., Garcia-Rojas, J., Peimbert, A., & Mesa-Delgado, 
A. 2009, ApJ, 700, 654 

Ferland, G. J. Korista, K.T. Verner, D.A. Ferguson, J.W. Kingdon, J.B. Verner, & E.M. 
1998, PASP, 110, 761 

Fields, B.D. & Sarkar, S. 2008, in C. Amsler et al, Phys. Lett. B, 667, 1 

Gonzalez Delgado, R. M., Cervino, M., Martins, L. P., Leitherer, C, & Hauschildt, P. H. 
2005, MNRAS, 357, 945 

Gruenwald, R., Steigman, G., & Viegas, S. M. 2002, ApJ, 567, 931 

Hogan, C.J., Olive, K.A., & Scully, S.T. 1997, ApJ, 489, L119 

Hosford, A., Ryan, S. G., Garcia Perez, A. E., Norris, J. E., & Olive, K. A. 2009, A&A, 493, 
601 

Hummer, D. G. & Storey, P. J. 1987, MNRAS, 224, 801 

Izotov, Y. I., Chaffee, F. H., Foltz, C. B., Green, R. F., Guseva, N. G., & Thuan, T. X. 1999, 
ApJ, 527, 757 

Izotov, Y. I., & Thuan, T. X. 1998, ApJ, 497, 227 

Izotov, Y. I., & Thuan, T. X. 1998, ApJ, 500, 188 

Izotov, Y. I., & Thuan, T. X. 2004, ApJ, 602, 200 (IT04) 

Izotov, Y. I., Thuan, T. X., & Lipovetsky, V. A. 1994, ApJ, 435, 647 

Izotov, Y. I., Thuan, T. X., & Lipovetsky, V. A. 1997, ApJS, 108, 1 

Izotov, Y. I., Thuan, T. X., & Stasinska, G. 2007, ApJ, 662, 15 (ITS07) 

Izotov, Y. I., & Thuan, T. X.. arXiv: 1001 .44401 [astro-ph.CO] 

Kirkman, D., Tytler, D., Suzuki, N., O'Meara, J. M., & Lubin, D. 2003, ApJS, 149, 1 
Komatsu, E., et al. 2009, ApJS, 180, 330 

Luridiana, V., Peimbert, A., Peimbert, M., & Cervino, M. 2003, ApJ, 592, 846 



Malaney, R. A. k Mathews, G. J. 1993, Phys. Rep., 229, 145 

Martins, L. P., Gonzalez Delgado, R. M., Leitherer, C, Cervino, M., k Hauschildt, P. 2005, 
MNRAS, 358, 49 

McQuinn, K. B. W., Skillman, E. D., Cannon, J. M., Dalcanton, J. J., Dolphin, A., Stark, 
D., k Weisz, D. 2009, ApJ, 695, 561 

Olive, K. A., k Skillman, E. D. 2001, New Astronomy, 6, 119 (Paper I) 

Olive, K. A., k Skillman, E. D. 2004, ApJ, 617, 29 (Paper II) 

Olive, K. A., Skillman, E. D., k Steigman, G. 1997, ApJ, 483, 788 

Olive, K. A., Steigman, G. k Walker, T.P 2000, Phys Rep, 333, 389 

Pagel, B. E. J., Simonson, E. A., Terlevich, R. J., & Edmunds, M. G. 1992, MNRAS, 255, 
325 

Peimbert, M. 1967, ApJ, 150, 825 

Peimbert, M., Luridiana, V., k Peimbert, A. 2007, ApJ, 666, 636 (PLP07) 
Peimbert, A., Peimbert, M., k Luridiana, V. 2002, ApJ, 565, 668 
Peimbert, M., Peimbert, A., k Ruiz, M. T. 2000, ApJ, 541, 688 
Peimbert, M., k Torres-Peimbert, S. 1974, ApJ, 193, 327 

Pettini, M., Zych, B. J., Murphy, M. T., Lewis, A., k Steidel, C. C. 2008, MNRAS, 391, 
1499 

Porter, R. L., Bauman, R. P., Ferland, G. J. , k MacAdam, K. B. 2005, ApJ, 622L, 73 
Porter, R. L., Ferland, G. J., k MacAdam, K. B. 2007, ApJ, 657, 327 (PFM) 
Porter, R. L., Ferland, G. J., MacAdam, K. B., k Storey, P. J. 2009, MNRAS, 393, L36 
Sarkar, S. 1996 Rep. Prog. Phys., 59, 1493 
Seaton, M. J. 1979, MNRAS, 187, 73p 



-51 - 

Serpico, P. D., Esposito, S., Iocco, F., Mangano, G., Miele, G., & Pisanti, O. 2004, Journal 
of Cosmology and Astro-Particle Physics, 12, 10 

Skillman, E. D. 1985, ApJ, 290, 449 

Skillman, E. D. & Kennicutt, R. C. 1993, ApJ, 411, 655 

Skillman, E. D., Terlevich, E., & Terlevich, R. 1998, Space Science Reviews, 84, 105 

Spergel, D. N., Verde, L., Peiris, H. V., et al. 2003, ApJS, 148, 175 

Spergel, D. N., Bean, R., Dore, O., et al. 2006, ApJS, 170, 377 

Stasihska, G., & Izotov, Y. I. 2001, A&A, 378, 817 

Steigman, G., Viegas, S. M., & Gruenwald, R. 1997, ApJ, 490, 187 

Thuan, T. X., & Izotov, Y. I. 2002, Space Sci. Reviews, 100, 263 

Viegas, S. M., Gruenwald, R., & Steigman, G. 2000, ApJ, 531, 813 

Walker, T.P, Steigman, G., Schramm, D.N., Olive, K.A., & Kang, K. 1991, ApJ, 376, 51 



This preprint was prepared with the A AS IATgX macros v5.2. 



-52 - 



Table 8. Primordial Helium Regression Values 



Object 


Hc+/H+ 


Hc++/H+ 
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O/H x 10 4 


I Zw 18 


0.08008 ± 0.00469 


0.0024 ± 0.0024 


0.2480 ± 0.0119 


0.15 ± 0.01 


SBS 0335-052 


0.08536 ± 0.00644 


0.0023 ± 0.0023 


0.2595 ± 0.0150 


0.19 ± 0.01 


SBS 0940+544N 


0.08506 ± 0.00484 


0.0000 


0.2537 ± 0.0108 


0.26 ± 0.01 


SBS 1159+545 


0.08332 ± 0.00390 


0.0010 ± 0.0010 


0.2520 ± 0.0090 


0.30 ± 0.01 


SBS 1420+544 


0.08939 ± 0.00393 


0.0011 ± 0.0011 


0.2655 ± 0.0088 


0.54 ± 0.01 


Haro 29 


0.08576 ± 0.00174 


0.0011 ± 0.0011 


0.2575 ± 0.0045 


0.58 ± 0.01 


Mrk 193 


0.08660 ± 0.00435 


0.0010 ± 0.0010 


0.2592 ± 0.0098 


0.61 ± 0.02 


NGC 2363A 


0.08328 ± 0.00244 


0.0012 ± 0.0012 


0.2522 ± 0.0061 


0.69 ± 0.01 


NGC 346 


0.08777 ± 0.00168 


0.0000 


0.2593 ± 0.0037 


1.02 ± 0.02 



